Greedy reduction of Bacillus subtilis genome yields emergent phenotypes of high resistance to a DNA damaging agent and low evolvability.

Genome-scale engineering enables rational removal of dispensable genes in chassis genomes. Deviating from this approach, we applied greedy accumulation of deletions of large dispensable regions in the Bacillus subtilis genome, yielding a library of 298 strains with genomes reduced up to 1.48 Mb in size. High-throughput physiological phenotyping of these strains confirmed that genome reduction is associated with substantial loss of cell fitness and accumulation of synthetic-sick interactions. Transcriptome analysis indicated that <15% of the genes conserved in our genome-reduced strains exhibited a twofold or higher differential expression and revealed a thiol-oxidative stress response. Most transcriptional changes can be explained by loss of known functions and by aberrant transcription at deletion boundaries. Genome-reduced strains exhibited striking new phenotypes relative to wild type, including a very high resistance (increased >300-fold) to the DNA-damaging agent mitomycin C and a very low spontaneous mutagenesis (reduced 100-fold). Adaptive laboratory evolution failed to restore cell fitness, except when coupled with a synthetic increase of the mutation rate, confirming low evolvability. Although mechanisms underlying this emergent phenotype are not understood, we propose that low evolvability can be leveraged in an engineering strategy coupling reductive cycles with evolutive cycles under induced mutagenesis.

Termination factor Rho mediates transcriptional reprogramming of Bacillus subtilis stationary phase.

Transcription termination factor Rho is known for its ubiquitous role in suppression of pervasive, mostly antisense, transcription. In the model Gram-positive bacterium Bacillus subtilis, de-repression of pervasive transcription by inactivation of rho revealed the role of Rho in the regulation of post-exponential differentiation programs. To identify other aspects of the regulatory role of Rho during adaptation to starvation, we have constructed a B. subtilis strain (Rho+) that expresses rho at a relatively stable high level in order to compensate for its decrease in the wild-type cells entering stationary phase. The RNAseq analysis of Rho+, WT and Δrho strains (expression profiles can be visualized at http://genoscapist.migale.inrae.fr/seb_rho/) shows that Rho over-production enhances the termination efficiency of Rho-sensitive terminators, thus reducing transcriptional read-through and antisense transcription genome-wide. Moreover, the Rho+ strain exhibits global alterations of sense transcription with the most significant changes observed for the AbrB, CodY, and stringent response regulons, forming the pathways governing the transition to stationary phase. Subsequent physiological analyses demonstrated that maintaining rho expression at a stable elevated level modifies stationary phase-specific physiology of B. subtilis cells, weakens stringent response, and thereby negatively affects the cellular adaptation to nutrient limitations and other stresses, and blocks the development of genetic competence and sporulation. These results highlight the Rho-specific termination of transcription as a novel element controlling stationary phase. The release of this control by decreasing Rho levels during the transition to stationary phase appears crucial for the functionality of complex gene networks ensuring B. subtilis survival in stationary phase.

Two functionally distinct heme/iron transport systems are virulence determinants of the fish pathogen Flavobacterium psychrophilum.

Bacterial pathogens have a critical impact on aquaculture, a sector that accounts for half of the human fish consumption. Flavobacterium psychrophilum (phylum Bacteroidetes) is responsible for bacterial cold-water disease in salmonids worldwide. The molecular factors involved in host invasion, colonization and haemorrhagic septicaemia are mostly unknown. In this study, we identified two new TonB-dependent receptors, HfpR and BfpR, that are required for adaptation to iron conditions encountered during infection and for virulence in rainbow trout. Transcriptional analyses revealed that their expression is tightly controlled and upregulated under specific iron sources and concentrations. Characterization of deletion mutants showed that they act without redundancy: BfpR is required for optimal growth in the presence of high haemoglobin level, while HfpR confers the capacity to acquire nutrient iron from haem or haemoglobin under iron scarcity. The gene hfpY, co-transcribed with hfpR, encodes a protein related to the HmuY family. We demonstrated that HfpY binds haem and contributes significantly to host colonization and disease severity. Overall, these results are consistent with a model in which both BfpR and Hfp systems promote haem uptake and respond to distinct signals to adapt iron acquisition to the different stages of pathogenesis. Our findings give insight into the molecular basis of pathogenicity of a serious pathogen belonging to the understudied family Flavobacteriaceae and point to the newly identified haem receptors as promising targets for antibacterial development.

New genetic biomarkers to differentiate non-pathogenic from clinically relevant Bacillus cereus strains.

OBJECTIVES: Bacillus cereus is responsible for food poisoning and rare but severe clinical infections. The pathogenicity of strains varies from harmless to lethal strains. However, there are currently no markers, either alone or in combination, to differentiate pathogenic from non-pathogenic strains. The objective of the study was to identify new genetic biomarkers to differentiate non-pathogenic from clinically relevant B. cereus strains. METHODS: A first set of 15 B. cereus strains were compared by RNAseq. A logistic regression model with lasso penalty was applied to define combination of genes whose expression was associated with strain pathogenicity. The identified markers were checked for their presence/absence in a collection of 95 B. cereus strains with varying pathogenic potential (food-borne outbreaks, clinical and non-pathogenic). Receiver operating characteristic area under the curve (AUC) analysis was used to determine the combination of biomarkers, which best differentiate between the "disease" versus "non-disease" groups. RESULTS: Seven genes were identified during the RNAseq analysis with a prediction to differentiate between pathogenic and non-pathogenic strains. The validation of the presence/absence of these genes in a larger collection of strains coupled with AUC prediction showed that a combination of four biomarkers was sufficient to accurately discern clinical strains from harmless strains, with an AUC of 0.955, sensitivity of 0.9 and specificity of 0.86. CONCLUSIONS: These new findings help in the understanding of B. cereus pathogenic potential and complexity and may provide tools for a better assessment of the risks associated with B. cereus contamination to improve patient health and food safety.

Implication of a Key Region of Six Bacillus cereus Genes Involved in Siroheme Synthesis, Nitrite Reductase Production and Iron Cluster Repair in the Bacterial Response to Nitric Oxide Stress.

Bacterial response to nitric oxide (NO) is of major importance for bacterial survival. NO stress is a main actor of the eukaryotic immune response and several pathogenic bacteria have developed means for detoxification and repair of the damages caused by NO. However, bacterial mechanisms of NO resistance by Gram-positive bacteria are poorly described. In the opportunistic foodborne pathogen Bacillus cereus, genome sequence analyses did not identify homologs to known NO reductases and transcriptional regulators, such as NsrR, which orchestrate the response to NO of other pathogenic or non-pathogenic bacteria. Using a transcriptomic approach, we investigated the adaptation of B. cereus to NO stress. A cluster of 6 genes was identified to be strongly up-regulated in the early phase of the response. This cluster contains an iron-sulfur cluster repair enzyme, a nitrite reductase and three enzymes involved in siroheme biosynthesis. The expression pattern and close genetic localization suggest a functional link between these genes, which may play a pivotal role in the resistance of B. cereus to NO stress during infection.

Genoscapist: online exploration of quantitative profiles along genomes via interactively customized graphical representations.

SUMMARY: Genoscapist is a tool to design web interfaces generating high-quality images for interactive visualization of hundreds of quantitative profiles along a reference genome together with various annotations. Relevance is demonstrated by deployment of two websites dedicated to large condition-dependent transcriptome datasets available for Bacillus subtilis and Staphylococcus aureus. AVAILABILITY AND IMPLEMENTATION: Websites and source code freely accessible at https://genoscapist.migale.inrae.fr.

Transcriptome architecture and regulation at environmental transitions in flavobacteria: the case of an important fish pathogen.

The family Flavobacteriaceae (phylum Bacteroidetes) is a major component of soil, marine and freshwater ecosystems. In this understudied family, Flavobacterium psychrophilum is a freshwater pathogen that infects salmonid fish worldwide, with critical environmental and economic impact. Here, we report an extensive transcriptome analysis that established the genome map of transcription start sites and transcribed regions, predicted alternative sigma factor regulons and regulatory RNAs, and documented gene expression profiles across 32 biological conditions mimicking the pathogen life cycle. The results link genes to environmental conditions and phenotypic traits and provide insights into gene regulation, highlighting similarities with better known bacteria and original characteristics linked to the phylogenetic position and the ecological niche of the bacterium. In particular, osmolarity appears as a signal for transition between free-living and within-host programs and expression patterns of secreted proteins shed light on probable virulence factors. Further investigations showed that a newly discovered sRNA widely conserved in the genus, Rfp18, is required for precise expression of proteases. By pointing proteins and regulatory elements probably involved in host-pathogen interactions, metabolic pathways, and molecular machineries, the results suggest many directions for future research; a website is made available to facilitate their use to fill knowledge gaps on flavobacteria.

Exploring Listeria monocytogenes Transcriptomes in Correlation with Divergence of Lineages and Virulence as Measured in Galleria mellonella.

As for many opportunistic pathogens, the virulence potential of Listeria monocytogenes is highly heterogeneous between isolates and correlated, to some extent, with phylogeny and gene repertoires. In sharp contrast with copious data on intraspecies genome diversity, little is known about transcriptome diversity despite the role of complex genetic regulation in pathogenicity. The current study implemented RNA sequencing to characterize the transcriptome profiles of 33 isolates under optimal in vitro growth conditions. Transcript levels of conserved single-copy genes were comprehensively explored from several perspectives, including phylogeny, in silico-predicted virulence category based on epidemiological multilocus sequence typing (MLST) data, and in vivo virulence phenotype assessed in Galleria mellonella Comparing baseline transcriptomes between isolates was intrinsically more complex than standard genome comparison because of the inherent plasticity of gene expression in response to environmental conditions. We show that the relevance of correlation analyses and their statistical power can be enhanced by using principal-component analysis to remove the first level of irrelevant, highly coordinated changes linked to growth phase. Our results highlight the major contribution of transcription factors with key roles in virulence to the diversity of transcriptomes. Divergence in the basal transcript levels of a substantial fraction of the transcriptome was observed between lineages I and II, echoing previously reported epidemiological differences. Correlation analysis with in vivo virulence identified numerous sugar metabolism-related genes, suggesting that specific pathways might play roles in the onset of infection in G. mellonellaIMPORTANCEListeria monocytogenes is a multifaceted bacterium able to proliferate in a wide range of environments from soil to mammalian host cells. The accumulated genomic data underscore the contribution of intraspecies variations in gene repertoire to differential adaptation strategies between strains, including infection and stress resistance. It seems very likely that the fine-tuning of the transcriptional regulatory network is also a key component of the phenotypic diversity, albeit more difficult to investigate than genome content. Some studies reported incongruity in the basal transcriptome between isolates, suggesting a putative relationship with phenotypes, but small isolate numbers hampered proper correlation analyses with respect to their characteristics. The present study is the embodiment of the promising approach that consists of analyzing correlations between transcriptomes and various isolate characteristics. Statistically significant correlations were found with phylogenetic groups, epidemiological evidence of virulence potential, and virulence in Galleria mellonella larvae used as an in vivo model.

Bacillus cereus, a serious cause of nosocomial infections: Epidemiologic and genetic survey.

Bacillus cereus is the 2nd most frequent bacterial agent responsible for food-borne outbreaks in France and the 3rd in Europe. In addition, local and systemic infections have been reported, mainly describing individual cases or single hospital setting. The real incidence of such infection is unknown and information on genetic and phenotypic characteristics of the incriminated strains is generally scarce. We performed an extensive study of B. cereus strains isolated from patients and hospital environments from nine hospitals during a 5-year study, giving an overview of the consequences, sources and pathogenic patterns of B. cereus clinical infections. We demonstrated the occurrence of several hospital-cross-contaminations. Identical B. cereus strains were recovered from different patients and hospital environments for up to 2 years. We also clearly revealed the occurrence of inter hospital contaminations by the same strain. These cases represent the first documented events of nosocomial epidemy by B. cereus responsible for intra and inter hospitals contaminations. Indeed, contamination of different patients with the same strain of B. cereus was so far never shown. In addition, we propose a scheme for the characterization of B. cereus based on biochemical properties and genetic identification and highlight that main genetic signatures may carry a high pathogenic potential. Moreover, the characterization of antibiotic resistance shows an acquired resistance phenotype for rifampicin. This may provide indication to adjust the antibiotic treatment and care of patients.

Genomic Diversity and Evolution of the Fish Pathogen Flavobacterium psychrophilum.

Flavobacterium psychrophilum, the etiological agent of rainbow trout fry syndrome and bacterial cold-water disease in salmonid fish, is currently one of the main bacterial pathogens hampering the productivity of salmonid farming worldwide. In this study, the genomic diversity of the F. psychrophilum species is analyzed using a set of 41 genomes, including 30 newly sequenced isolates. These were selected on the basis of available MLST data with the two-fold objective of maximizing the coverage of the species diversity and of allowing a focus on the main clonal complex (CC-ST10) infecting farmed rainbow trout (Oncorhynchus mykiss) worldwide. The results reveal a bacterial species harboring a limited genomic diversity both in terms of nucleotide diversity, with ~0.3% nucleotide divergence inside CDSs in pairwise genome comparisons, and in terms of gene repertoire, with the core genome accounting for ~80% of the genes in each genome. The pan-genome seems nevertheless "open" according to the scaling exponent of a power-law fitted on the rate of new gene discovery when genomes are added one-by-one. Recombination is a key component of the evolutionary process of the species as seen in the high level of apparent homoplasy in the core genome. Using a Hidden Markov Model to delineate recombination tracts in pairs of closely related genomes, the average recombination tract length was estimated to ~4.0 Kbp and the typical ratio of the contributions of recombination and mutations to nucleotide-level differentiation (r/m) was estimated to ~13. Within CC-ST10, evolutionary distances computed on non-recombined regions and comparisons between 22 isolates sampled up to 27 years apart suggest a most recent common ancestor in the second half of the nineteenth century in North America with subsequent diversification and transmission of this clonal complex coinciding with the worldwide expansion of rainbow trout farming. With the goal to promote the development of tools for the genetic manipulation of F. psychrophilum, a particular attention was also paid to plasmids. Their extraction and sequencing to completion revealed plasmid diversity that remained hidden to classical plasmid profiling due to size similarities.

PamR, a new MarR-like regulator affecting prophages and metabolic genes expression in Bacillus subtilis.

B. subtilis adapts to changing environments by reprogramming its genetic expression through a variety of transcriptional regulators from the global transition state regulators that allow a complete resetting of the cell genetic expression, to stress specific regulators controlling only a limited number of key genes required for optimal adaptation. Among them, MarR-type transcriptional regulators are known to respond to a variety of stresses including antibiotics or oxidative stress, and to control catabolic or virulence gene expression. Here we report the characterization of the ydcFGH operon of B. subtilis, containing a putative MarR-type transcriptional regulator. Using a combination of molecular genetics and high-throughput approaches, we show that this regulator, renamed PamR, controls directly its own expression and influence the expression of large sets of prophage-related and metabolic genes. The extent of the regulon impacted by PamR suggests that this regulator reprograms the metabolic landscape of B. subtilis in response to a yet unknown signal.

Termination factor Rho: From the control of pervasive transcription to cell fate determination in Bacillus subtilis.

In eukaryotes, RNA species originating from pervasive transcription are regulators of various cellular processes, from the expression of individual genes to the control of cellular development and oncogenesis. In prokaryotes, the function of pervasive transcription and its output on cell physiology is still unknown. Most bacteria possess termination factor Rho, which represses pervasive, mostly antisense, transcription. Here, we investigate the biological significance of Rho-controlled transcription in the Gram-positive model bacterium Bacillus subtilis. Rho inactivation strongly affected gene expression in B. subtilis, as assessed by transcriptome and proteome analysis of a rho-null mutant during exponential growth in rich medium. Subsequent physiological analyses demonstrated that a considerable part of Rho-controlled transcription is connected to balanced regulation of three mutually exclusive differentiation programs: cell motility, biofilm formation, and sporulation. In the absence of Rho, several up-regulated sense and antisense transcripts affect key structural and regulatory elements of these differentiation programs, thereby suppressing motility and biofilm formation and stimulating sporulation. We dissected how Rho is involved in the activity of the cell fate decision-making network, centered on the master regulator Spo0A. We also revealed a novel regulatory mechanism of Spo0A activation through Rho-dependent intragenic transcription termination of the protein kinase kinB gene. Altogether, our findings indicate that distinct Rho-controlled transcripts are functional and constitute a previously unknown built-in module for the control of cell differentiation in B. subtilis. In a broader context, our results highlight the recruitment of the termination factor Rho, for which the conserved biological role is probably to repress pervasive transcription, in highly integrated, bacterium-specific, regulatory networks.

Large-scale reduction of the Bacillus subtilis genome: consequences for the transcriptional network, resource allocation, and metabolism.

Understanding cellular life requires a comprehensive knowledge of the essential cellular functions, the components involved, and their interactions. Minimized genomes are an important tool to gain this knowledge. We have constructed strains of the model bacterium, Bacillus subtilis, whose genomes have been reduced by ∼36%. These strains are fully viable, and their growth rates in complex medium are comparable to those of wild type strains. An in-depth multi-omics analysis of the genome reduced strains revealed how the deletions affect the transcription regulatory network of the cell, translation resource allocation, and metabolism. A comparison of gene counts and resource allocation demonstrates drastic differences in the two parameters, with 50% of the genes using as little as 10% of translation capacity, whereas the 6% essential genes require 57% of the translation resources. Taken together, the results are a valuable resource on gene dispensability in B. subtilis, and they suggest the roads to further genome reduction to approach the final aim of a minimal cell in which all functions are understood.

Staphylococcus aureus Transcriptome Architecture: From Laboratory to Infection-Mimicking Conditions.

Staphylococcus aureus is a major pathogen that colonizes about 20% of the human population. Intriguingly, this Gram-positive bacterium can survive and thrive under a wide range of different conditions, both inside and outside the human body. Here, we investigated the transcriptional adaptation of S. aureus HG001, a derivative of strain NCTC 8325, across experimental conditions ranging from optimal growth in vitro to intracellular growth in host cells. These data establish an extensive repertoire of transcription units and non-coding RNAs, a classification of 1412 promoters according to their dependence on the RNA polymerase sigma factors SigA or SigB, and allow identification of new potential targets for several known transcription factors. In particular, this study revealed a relatively low abundance of antisense RNAs in S. aureus, where they overlap only 6% of the coding genes, and only 19 antisense RNAs not co-transcribed with other genes were found. Promoter analysis and comparison with Bacillus subtilis links the small number of antisense RNAs to a less profound impact of alternative sigma factors in S. aureus. Furthermore, we revealed that Rho-dependent transcription termination suppresses pervasive antisense transcription, presumably originating from abundant spurious transcription initiation in this A+T-rich genome, which would otherwise affect expression of the overlapped genes. In summary, our study provides genome-wide information on transcriptional regulation and non-coding RNAs in S. aureus as well as new insights into the biological function of Rho and the implications of spurious transcription in bacteria.

Constitutive Stringent Response Restores Viability of Bacillus subtilis Lacking Structural Maintenance of Chromosome Protein.

Bacillus subtilis mutants lacking the SMC-ScpAB complex are severely impaired for chromosome condensation and partitioning, DNA repair, and cells are not viable under standard laboratory conditions. We isolated suppressor mutations that restored the capacity of a smc deletion mutant (Δsmc) to grow under standard conditions. These suppressor mutations reduced chromosome segregation defects and abrogated hypersensitivity to gyrase inhibitors of Δsmc. Three suppressor mutations were mapped in genes involved in tRNA aminoacylation and maturation pathways. A transcriptomic survey of isolated suppressor mutations pointed to a potential link between suppression of Δsmc and induction of the stringent response. This link was confirmed by (p)ppGpp quantification which indicated a constitutive induction of the stringent response in multiple suppressor strains. Furthermore, sublethal concentrations of arginine hydroxamate (RHX), a potent inducer of stringent response, restored growth of Δsmc under non permissive conditions. We showed that production of (p)ppGpp alone was sufficient to suppress the thermosensitivity exhibited by the Δsmc mutant. Our findings shed new light on the coordination between chromosome dynamics mediated by SMC-ScpAB and other cellular processes during rapid bacterial growth.

Deciphering Genome Content and Evolutionary Relationships of Isolates from the Fungus Magnaporthe oryzae Attacking Different Host Plants.

Deciphering the genetic bases of pathogen adaptation to its host is a key question in ecology and evolution. To understand how the fungus Magnaporthe oryzae adapts to different plants, we sequenced eight M. oryzae isolates differing in host specificity (rice, foxtail millet, wheat, and goosegrass), and one Magnaporthe grisea isolate specific of crabgrass. Analysis of Magnaporthe genomes revealed small variation in genome sizes (39-43 Mb) and gene content (12,283-14,781 genes) between isolates. The whole set of Magnaporthe genes comprised 14,966 shared families, 63% of which included genes present in all the nine M. oryzae genomes. The evolutionary relationships among Magnaporthe isolates were inferred using 6,878 single-copy orthologs. The resulting genealogy was mostly bifurcating among the different host-specific lineages, but was reticulate inside the rice lineage. We detected traces of introgression from a nonrice genome in the rice reference 70-15 genome. Among M. oryzae isolates and host-specific lineages, the genome composition in terms of frequencies of genes putatively involved in pathogenicity (effectors, secondary metabolism, cazome) was conserved. However, 529 shared families were found only in nonrice lineages, whereas the rice lineage possessed 86 specific families absent from the nonrice genomes. Our results confirmed that the host specificity of M. oryzae isolates was associated with a divergence between lineages without major gene flow and that, despite the strong conservation of gene families between lineages, adaptation to different hosts, especially to rice, was associated with the presence of a small number of specific gene families. All information was gathered in a public database (http://genome.jouy.inra.fr/gemo).

Deciphering the molecular basis of wine yeast fermentation traits using a combined genetic and genomic approach.

The genetic basis of the phenotypic diversity of yeast is still poorly understood. Wine yeast strains have specific abilities to grow and ferment under stressful conditions compared with other strains, but the genetic basis underlying these traits is unknown. Understanding how sequence variation influences such phenotypes is a major challenge to address adaptation mechanisms of wine yeast. We aimed to identify the genetic basis of fermentation traits and gain insight into their relationships with variations in gene expression among yeast strains. We combined fermentation trait QTL mapping and expression profiling of fermenting cells in a segregating population from a cross between a wine yeast derivative and a laboratory strain. We report the identification of QTL for various fermentation traits (fermentation rates, nitrogen utilization, metabolites production) as well as expression QTL (eQTL). We found that many transcripts mapped to several eQTL hotspots and that two of them overlapped with QTL for fermentation traits. A QTL controlling the maximal fermentation rate and nitrogen utilization overlapping with an eQTL hotspot was dissected. We functionally demonstrated that an allele of the ABZ1 gene, localized in the hotspot and involved in p-aminobenzoate biosynthesis, controls the fermentation rate through modulation of nitrogen utilization. Our data suggest that the laboratory strain harbors a defective ABZ1 allele, which triggers strong metabolic and physiological alterations responsible for the generation of the eQTL hotspot. They also suggest that a number of gene expression differences result from some alleles that trigger major physiological disturbances.